Hexagonally Warped Dirac Cones and Topological Phase Transition 

in Silicene Superstructure 
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Silicene is a monolayer of silicon atoms forming a two-dimensional honeycomb lattice. We investigate elec- 
tronic properties of the (4x4) superstructure, which is formed on top of the Ag(lll) substrate. The band 
structure is obtained by folding that of silicene according to the superstructure. It is a peculiar feature that the 
two inequivalent K and K' points in the silicene Brillouin zone are identified by the folding. Consequently, 
two Dirac cones coexist at the same Dirac point in the momentum space. They are hexagonally warped by the 
Coulomb interaction. The silicene superstructure is a quantum spin-Hall insulator. It is surprising that the band 
gap never closes even if we apply the electric field. It is highly contrasted to the case of free-standing silicene, 
where a topological phase transition occurs. We then study generic superstructure and derive the conditions 
when the K and K' points are identified and when a topological phase transition may occur. 



Graphene has attracted tremendous attention. Silicene is a 
sheet of silicon atoms replacing carbons in graphene. It could 
follow this trend llHlJ] and open new perspectives for appli- 
cations, especially due to its compatibility with Si-based elec- 
tronics. Almost every striking property of graphene could be 
transferred to this innovative material. Indeed, it has Dirac 
cones akin to graphene. It has additionally a salient feature, 
that is a relatively large spin-orbit (SO) gap[6], which pro- 
vides a mass to Dirac electrons and realizes a detectable quan- 
tum spin Hall (QSH) effect0@J. Furthermore a topological 
phase transition occurs from the QSH insulator to the trivial 
band insulator in the electric field|9]. However, these prop- 
erties are those of free-standing silicene. Silicene has been 
synthesized on metal surfaces that interact weakly with the 
silicon atoms llHSI]. 

Silicene available now is mostly grown on the Ag(lll) 
substrate. The most often observed one is the so-called 
(4x4) superstructure |H Q], though many other superstruc- 
tures are possible due to the commensurability with the Ag 
substrate! 13, Note that the symbol (4x4) refers to the 
structure of the Ag substrate. Indeed, there are 3 x 3 sili- 
con atoms on top of 4 x 4 silver atoms in this structure. We 
should call it the (4x4) Ag- superstructure or the (3x3) Si- 
superstructure. The unit cell contains 2x(3x3) = 18 silicon 
atoms. The superstructure contain two sublattices: There are 
6 (12) atoms in the higher (lower) sublattice which is 3. OA 
(2. 2 A) above the Ag surface, as illustrated in FigQ] This 
structure has been observed in scanning tunneling microscopy 
(STM) images, where only the atoms in the higher sublattice 
are visibleflll. 

In silicene the states near the Fermi energy are tt orbitals 
residing near the inequivalent K and K' points at opposite cor- 
ners of the hexagonal Brillouin zone. We refer to them also 
as the point with the valley index r] = ±1. It is folded 
into a reduced Brillouin zone in a superstructure. Hence we 
can generally argue that the low-energy theory is described by 
massive Dirac electrons and that a silicene superstructure is 
a QSH insulator. However it is a nontrivial problem whether 
there exists a topological phase transition in the external elec- 
tric field(9|]. Furthermore, we may expect an intriguing phe- 




FIG. 1 : (Color online) Unit cell of the (3 x 3) Si-superstructure. There 
are 18 silicon atoms per unit cell, in which 6 atoms are in higher 
position and 12 atoms are in lower position. There are 3 A and 3 B 
sites in higher position. Their vertical sepation is 2£ « 0.8 A.. The 
basis vectors are represented by a^. 



nomenon where the K and K' points are identified by the fold- 
ing. When it happens, there exists only one Dirac point where 
two Dirac cones coexist, and this is experimentally detectable. 

In this paper we investigate electronic properties of the 
(3x3) Si- superstructure in detail. We show that the K and 
K' points are identified and that the band structure exhibits 
hexagonal warping in the presence of the Coulomb interac- 
tion. We also show that the band gap never closes even if we 
apply the electric field. We also make a discussion on generic 
silicene superstructure if the K and K' points are identified 
and if a topological phase transition occurs. 

The (3x3) Si- superstructure system is described by the 
36-band second-nearest-neighbor tight-binding model con- 
structed as follows. The basic term is[ 15] 



(i,j)a 
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where c\ a creates an electron with spin polarization a at 
site i, and / run over all the nearest/next-nearest 

neighbor hopping sites. The first term represents the usual 
nearest-neighbor hopping with the transfer energy t = 1.6eV. 
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FIG. 2: (Color online) Brillouin zone of the (3x3) Si-superstructure. 
Its shape is rhombus. It is constructed by folding the Brillouin zone 
of silicene. The K and K' points are identified, which we denote by 
K and take at the center of the Brillouin zone. The band structure 
along the horizontal and vertical thin lines are shown in Fig(3] 

The second term represents the intrinsic SO coupling with 
Aso = 3.9meV, where a = (cr Xj a yj a z ) is the Pauli ma- 
trix of spin, vij = (di x dj) / \di x dj\ with di and dj the 
two bonds connecting the next-nearest neighbors. Addition- 
ally there exist two types of Rashba SO couplings (Ari and 
Ar2) to describe siliceneLlO]. They are quite small, and play 
no important roles in the present analysis, as far as we have 
numerically checked. To avoid unnecessary complications, 
we neglect them. 

The key feature of the (3 x 3) Si- superstructure is that it is 
made of higher and lower sublattices, which contain 6 and 
12 silicon atoms per unit cell [Fig HI. Each silicon atom is 
weakly bound to the silver substrate. The effect is summarized 
as a chemical potential difference U between silicon atoms 
in higher and lower sublattices. Additionally we may apply 
the electric field E z perpendicular to the plane. It generates 
a staggered sublattice potential oc 2£E Z between them. The 
effective Hamiltonian is 

H 1 = (U + 2£E Z ) C K (2) 

iGhigh 

where the summation is taken over higher sites: U is a con- 
stant of the order of O.leV, and 2£ « 0.8 A according to the 
first-principle calculation tl3|]. 

The total Hamiltonian is H = Hq + Hi, where the spin s z 
is a good quantum number. Hence, the Hamiltonian is decom- 
posed into the spin sectors indexed by s z = ±1. It is enough 
to analyze the 18-band tight-binding model. We set E z = 
unless otherwise stated. 

The Brillouin zone of a superstructure is constructed from 
that of silicene by identifying two momentum vectors k and 
k' provided they are different only by the principal reciprocal 
vectors G\ and G 2 , as illustrated in FigEl See (IT2l) for the 
explicit form of Gi . The band structure satisfies 

e(k + Gi)=e(k). (3) 

Accordingly the band structure is simply constructed by fold- 
ing that of silicene when U = 0. 



(a1) < > u=0 (b1) < ► U=0.5t 




FIG. 3: (Color online) Band structure of the (3x3) Si- superstructure 
in extended Brillouin zone (al,bl) along the k x -axis and (a2,b2) 
along the k y -axis. The vertical axis is the energy in unit of t. The 
horizontal axis is momentum k. (a3,b3) Contour plot of two trigo- 
nally warped Dirac cones. We have taken U = for (al,a2,a3) and 
U = 0.5* for (bl,b2,b3). We have set A S0 = 0.2t for illustration. 

There may occur an intriguing phenomenon. In general, 
when there are two integers rti such that 

K-K' = n 1 G 1 +n 2 G 2 , (4) 

the K and K' points are identified. We call it the K point, and 
choose K = (0, 0). Two Dirac cones coexist at the K point. 

This happens indeed in the (3x3) Si-superstructure, where 
the Brillouin zone is obtained by folding that of silicene two 
times along the /c^-axis and three times along the /c y -axis. 
The area of the Brillouin zone is 1 /9 of that of silicene, and 
its shape is rhombus. We show numerically calculated band 
structures of the superstructure for U = and U ^ in FigO 
where we also show the contour plot of the band structures. 
There are two trigonally warped Dirac cones at the center of 
the Brillouin zone [FigOa3,b3)]. 

These two Dirac cones cross each other. The level cross- 
ing will be transformed into the level anticrossing when Dirac 
electrons interact among themselves. We show the contour 
plots of the Dirac cones calculated based on the tight-binding 
model in FiglH where we have introduced an infinitesimal 
on-site Coulomb interaction. The two Dirac cones are now 
hexagonally warped, as must be experimentally detectable. 

To understand this phenomenon analytically, we analyze 
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FIG. 4: (Color online) Contour plots of two hexagonally warped 
Dirac cones calculated based on the tight-binding model. They 
emerge when two crossing levels are transformed into two anticross- 
ing levels by an infinitesimal on-site Coulomb interaction in trigo- 
nally warped Dirac dones in Fig[3ta3). We have taken U = for 
simplicity. The Dirac cones are modified only slightly when U ^ 0. 



the low energy Dirac theory of the Hamiltonian (Q]), 



Hrj = hv ¥ (k x r x + rjkyTy) 



37T7]i/2 

f/ac / 2 



(5) 



where r = (r Xj r yj r z ) is the Pauli matrix of pseudospin as- 
sociated with the A and B sites. Two Dirac cones are labelled 
by rj. The way of trigonal warping is represented by the factor 
e 37rW 2 w ith the phase being opposite (rj — ±1) between the 
two types of electrons ll 1611 . Note that we cannot incorporate 
the potential term © into this simple Dirac theory. Hence we 
assume U = in the analytical treatment. 

When we introduce the on-site Coulomb interaction, the 
Hamiltonian and are mixed into 



H 



H21 #K' 



(6) 



where we have chosen the basis {c^, , c\ , }, and 
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with V representing the Coulomb energy. Up to the first order 
in V, the two bands are explicitly given by 

e\ = (hv F k) 2 + A^ + a 2 ± V2ahv F k Vl + cos 6(9, (10) 

where is the azimuthal angle in the momentum space, k x = 
k cos 0. We can check that these two band structures repro- 
duce the hexagonally warped Dirac cones in Fig 0] quite well. 




FIG. 5: (Color online) Band structure of nanorbibbon made of sil- 
icene superstructure. The vertical axis is the energy in unit of t. The 
horizontal axis is the momentum k. Two points k = and k = 2tt 
are identified, and gives the K point, (a) U = 0. (b) U = 0.5t. We 
have set Aso = 0.2t for illustration. Two lines connect the K point in 
the conduction and valence bands after winding once a circle of the 
whole Brillouin zone. 
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FIG. 6: (Color online) Band gap A as a function of the electric field 
E' z = E z + U/2L We show A for the (3 x 3) Si- superstructure in the 
solid curves and for free-standing silicene in the dashed lines. The 
band gap closes at E' CT = Xso/i + U for the free-standing silicene. 
The vertical axis is the energy in unit of t. We have set Aso = 0.2t 
for illustration. 



The silicene superstructure is an insulator. In order to show 
that it is a QSH insulator, we have numerically calculated the 
band structure of a nanoribbon geometry in Fig El Two zero- 
energy edge modes are clearly seen: Each line represents the 
two-fold degenerate helical modes propagating one edge. It 
implies that the silicene superstructure is a QSH insulator. The 
band structure of a nanoribbon with U = is just the one ob- 
tained by folding that of a silicene nanoribbon as in FigOa). It 
is modified for U 7^ 0, but the modification is only slight even 
when we take a quite large value for U, as shown in FigJ3b). 
The essential feature is that the edge state connects the con- 
duction and valence bands at the same K point. Nevertheless, 
the edge state is topologically protected since it winds a circle 
of the whole Brillouin zone as in FiglSJ 

We now switch on the external electric field E z in (|2]), and 
discuss the silicene superstructure in its presence. We show 
the band gap as a function of E z in Fig(6j The band gap is 
found to be almost insensitive to E z . It is surprising that the 
band gap never closes. This is highly contrasted to the case of 
free-standing silicene, where the band gap closes and a topo- 
logical phase transition occurs [9] at a certain critical field E CY . 
See (fTTl) for discussions on this point. 

We have studied the (3x3) Si- superstructure in detail. We 
proceed to make a discussion on generic silicene superstruc- 
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ture. All possible silicene superstructures are classified by the 
translational vectors g x and g 2 of the unit cell defined by 

g 1 =pa 1 + qa 2 , g 2 = e ±i7r/3 £i, (11) 

with the use of the two basis vectors a\ and a 2 of the honey- 
comb lattice [Fig[U, and < q < p. The choice of the angle 
±7r/3 generate a superstructure of type-I and type-II, respec- 
tively. When q = 0, type-I and type-II become identical. The 
reciprocal vectors G\ and G 2 are 

G i = (g i xn)/\g 1 xg 2 \, (12) 

where n is the unit vector perpendicular to the silicene plane. 
By a geometrical analysis of the honeycomb lattice we have 

\9i\ = Vp 2 + q 2 +pq- (13) 

The structure thus generated is the (|^ 1 | x |gf 1 |) Si- 
superstructure. The number of silicon atoms per unit cell is 
given by 

N = 2\ 9l \ 2 =2(p 2 + q 2 +pq). (14) 

The angle of the unit cell from the x-axis is given by 

tan0 = V3q/(2p + q). (15) 

It is another problem on which Ag- superstructure this Si- 
superstructure is grown. Here we give some correspondences 
between themlll3[ 1 1 411 : 



Ag- superstructure 


Si- superstructure 


P 


q 


ratio 


4 x 4[1, 2] 


3x3 


3 





1.33 


V7 xV7 


2x2 


2 





1.32 


VT3 x Vl3\2] 


V7 x V7 


2 


1 


1.36 


V21x V21 


2V3 x 2V3 


2 


2 


1.32 


2V3 x 2>/3[14] 


V7 x V7 


2 


1 


1.31 



(16) 



The ratios of the Ag and Si superstructures is around 1.3. 

We have found that there exists no topological phase transi- 
tion in the (3x3) Si- superstructure. It is important to explore 
the same problem in generic superstructure. We assume that 
the buckled structure consists of two sublattices. Indeed, there 
is a report! 13] that this is the case for all structures grown on 
Ag( 1,1,1). Let the numbers of A sites and B sites be Na and 
Nb, respectively, in the higher sublattice per unit cell which 
contains N silicon atoms. As we have stated, we cannot in- 
corporate the potential term into the simple Dirac theory 
(0). Nevertheless, it can be argued to yield a contribution to 
the Dirac mass m d within the mean-field approximation, 

Na - Np> 

m D = A N B (U + 2£E Z ) + t^Aso- (17) 

The phase transition would occur Hill at = 0, which 
determines the critical field E cv . In the case of free-standing 
silicene, we have U = 0, N A = 1, N B = 0, N = 2, and 
hence 

E cr = Xso/l, (18) 



as agrees with the previous result! 9]. In the present (3 x 3) Si- 
superstructure, we have Na = Nb = 3, N = 18, and we find 
E CY — » 00, or equivalently there is no phase transition point. 



In conclusion, we have investigated the (3x3) buckled su- 
perstructure of silicene. It is found that the K and K' points 
become identical and hence two Dirac cones coexist at the K 
point. They exhibit hexagonal warping in the presence of the 
Coulomb interaction, which must be observed by angle re- 
solved photoemission spectroscopy (ARPES). It is also found 
that there is no topological phase transition. Furthermore we 
have explored electronic properties of generic superstructure. 
The K and K' points become identified when the condition 
© is satisfied. We may expect a topological phase transition 
when the numbers of the A and B sites are different in the 
higher sublattice of the backed silicene superstructure. 
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